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Abstract In the present letter we attempt to find a new class of solutions relating 
radiation model by using Homotopy Perturbation Method (HPM). We consider 
that the spherical distribution is made of perfect fluid where radiation is along the 
radially outward direction of the system. Here we have employed HPM as a new 
tool for astrophysical systems which via the mass polynomial facilitates solving of 
field equations. A set of interior solutions has been found out on the basis of the 
equation of state in the form p = ip with the demand that the radiation model is 
consistent with the Schwarzschild metric. It is observed that the set of solutions 
thus developed provides new metric to explain strange star which however seems 
not free from singularity. 
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1 Introduction 

Under the general relativistic astrophysics the composition, nature and physical 
properties of compact stars have recently become a more fascinating research arena 
with diverse theoretical as well as applied issues. Several scientists fn[2i[3i have 
been pointed out that the quark matter, made of u, d and s quarks, which is 
known as the strange quark matter may be more stable than ordinary nuclear 
matter. This strange matter hypothesis actually indicates important properties 
of highly compact objects like neutron star and dwarf star. Essentially strange 
(quark) matter is a kind of Quark-Gluon-Plasma (QGP) which may be composed 
of equal numbers of u, d and s quarks along with a small number of electrons and 
some baryons. The strange stars (SS) named after the strange quark are therefore 
composed of strange matter. 

Using the MIT bag model Farook et al. [3] have obtained a deterministic model 
of strange star where they considered a mass polynomial and analyzed all physical 
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properties of strange star. However they were unable to show any properties up 
to 6 km from the centre of the system. 

In this present letter we have constructed a model to discuss the physical prop¬ 
erties of strange star altogether from the centre to the surface of the star. Solving 
the non-linear equations analytically has always been a challenge in astrophysics. 
The Homotopy Perturbation Method (HPM) is a powerful and very simple tool 
to solve these kind of equations with least number of assumptions. In the present 
work modeling the strange star we have developed expression of mass (mass poly¬ 
nomial) which is radial function of r. However, it is not just assumed arbitrarily 
rather we have done this by the help of HPM. Later on we have substituted that 
expression of mass to solve Einstein’s Field Equations. 

Therefore, following Rahaman et al. we have obtained the interior solution 
of the Einstein field equation of a spherically symmetric system of radiating star 
by introducing a little bit correction. In that paper [5] it was concluded that 
the model can not be compatible with the compact stars. In the present work 
however, with some appropriate changes, we have been able to justify the physical 
properties of the compact object, especially strange stars. We also would like to 
refer that the work of Rahaman et al. [5] where Homotopy Perturbation Method 
has been employed for a spherically symmetric system of radiating star suffers 
from instability problem. In the present letter our motivation is to present a stable 
model of a radiating compact star. 

The outline of our investigation is as follows: In the Sec. 2 we have calculated 
mass of the system using firstly, the Maximum Entropy Principle, and secondly, 
the Homotopy Perturbation Method. Sec. 3 deals with the solution of Einstein’s 
field equations for different physical parameters, viz. the pressures {pr and pt) and 
density. We have discussed and explored several physical features in Sec. 4 and a 
comparative study has been conducted in Sec. 5 for validity of the data set of the 
present model with that of the exiting strange stars available in the literature. In 
the last sec. 6 we have drawn some concluding remarks with some salient features 
of the present model. 


2 Calculation of mass of the spherical system 

2.1 The Maximum Entropy Principle 


The thermal radiation of a perfect fluid sphere satisfies following the equation of 
state (EOS) 


Pr = ip. 


( 1 ) 


It is to note that since the fluid is anisotropic so pt and pr are different but we 
have used Eq. © as our concern is the radiation model. 

The other necessary equations for the perfect fluid system are as follows: 



( 2 ) 



( 3 ) 
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Here s is the entropy density of the system, T is the temperature of the black 
body radiation, p is the matter density of the system and b is an arbitrary constant 
of order unity (in natural units G = c = h = k = 1). 

We consider the interior space-time metric of the spherical symmetric system 


as 


ds^ 


-9u{r)dl^ -h 


^ 2m(r) 

r 


dr^ + r^(d6'^ -|- sin^ 


(4) 


where m{r) is the mass distribntion of the system. The time-time component of 
the metric gu is the function of radial component r only. 

The general energy-momentum tensor for the spherically symmetric radiating 
perfect fluid system is as follows 


Tif = {p+ -I- (pt - ip)??'"??!/, 


(5) 


with u^Uu = —rj^Pv = 1. 

Here the vector is the fluid 4-velocity of the local rest frame and is the 
unit space-like vector of the radiating system,and the vectors are orthogonal to 
each other i.e. u^pi, = 0. 

The constraint relation comes from the time component of Einstein field equa¬ 
tion = S-TrT^i/ is 


m'[r) 

^ 47rj-2 


( 6 ) 


From Eq. ([2]), Eq. ([3|) and Eq. ® we have 


s = 





3 

4 


(7) 


Now one can write the total entropy of the spherical system for the matter 
distribution up to r < i? as 


S = 


s(r) 1 


2m(r) ^ 


- 1/2 

dV = a I Ldr, 
Jo 


( 8 ) 


where a = |47rfei and the Lagrangian of the system as 


/ /X 2 1 / 

L = (m ) ■* r 2 I 1 — 


2m(r) 


From the Euler-Lagrangian equation of motion 


(9) 



we therefore obtain 


3 // , 3 „ , 3 , 3 m'm 1 / /^2 ^ 

“ 16”" + S”" + S”" " 2 ^ ) = 0- 


( 10 ) 


( 11 ) 


The above Lagrangian equation of motion is a non-linear differential equa¬ 
tion. By solving this equation one would get the structure of the self gravitating 
radiation system. 
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2.2 The Homotopy Perturbation Method 

We have calculated the value of m by using Homotopy Perturbation Method (HPM)as 
provided by He [6] . 

The Homotopy for any non-linear differential equation is 

L{m) - L{mo)-\-p[L{mo) + N{m)\= 0, (12) 

where p is a parameter, such that p = 0 gives the initial approximation solution 
and p = 1 gives the desired solution for the mass. The term L(m) = ^ = m' is 
considered as the linear part of the equation and rest of the terms are non-linear. 
Also, it is adopted that L{mo) = mg. We choose initial solution for this system 
as mo = ar^, where a is a numerical constant. As a boundary condition we have 
assumed that the density should be minimum at the surface of the compact star so 
that d{p)/dr = 0 at the surface. One can see later on that after plotting different 
values of a (among three different values of a for a specific set of data), for a 
specihc value of a only, the density becomes minimum at the surface. 

To hnd out the value of m, we consider the general solution of m as 

m = mo-hp^mi-|-p^m 2 -h .... (13) 

Now from the Eq. (1111) . Eg. (1121) and Eq. (1131) we obtain 

m = limp^i{mo +p^mi + m 2 -|- ....) 

3 , 36 2 5 , 312 3 7 , n 2 , /I -x 

= ar +r -b-:r^o r -\-2>ac\r + C 2 (14) 

O oO 

where ci and C 2 are the integrating constants which will be evaluated later on in 
the next Sec. 3. 


3 The Solution of Einstein’s field equations 

The Einstein Eield Equation for the metric o, for the matter distribution given 
in Eq. m can be written as 



— - (l- ikl 

r3 V r ) gttr 



(15) 

(16) 




1 (g'u 

2 gtt 


gtt 


+ 


‘2r gtt 




m m'\ f I Ig'tt 


(17) 


From Eq. (USD, after substitution of Eq. (HI, we get the density of the system 


as 


P = 


47[-j.2 


477 


o I Of, 2 2 I 312 3 4 6aci 

3a -|- 3ba r H- —a r - 

5 r 


(18) 
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Fig. 1 Variation of mass as a function of radial distance r for a = —0.1233 x 10 ^ km ^ and 
R = 9.13 km 



Fig. 2 Variation of density as a function of radial distance r for a = —0.1233 x 10 ^ km ^ 
and R = 9.13 km 


Let US now choose p' = 0 at the surface (r = R) of the star, i.e. the density 
is minimum at the surface. By imposing this boundary condition in Eq. m , we 
obtain 

ons ^ _ 

(19) 


19 4- 2 p5 

Cl = IZaH H- —a H , 

5 


and from the boundary condition m = 0 at the centre r = 0, we get 

C2 = 0. (20) 

Therefore the mass of the spherical system, due to Eqs. (HI, (Hi) and dai, 
becomes 


m = ar'^ H—+ ?>ar {\2alt H- —aR°). (21) 

0 oO 0 


Hence Eq. Hi, by the help of Hi, takes the form 
1 


P = 


4 ?! 


3a + mar^ + ^ar + — (12ai?'^ + -^aR^) 
or 5 


( 22 ) 


The behaviour of the mass and density are shown in Eigs. 1 and 2. 
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Now from the EOS ([T]) and Eq. (I22II one get 




, ID 2 2 I 3 4 I d3 I 2o5n 

a-\-12a r H ——a r ^-(Iza/t H——a it ) 

5 r 5 


(23) 


5tt = O- 




(i-^y 


(24) 


where {r) 
dition, Eq. 


5 / d,r. After evaluating C, based on suitable boundary con- 

can be written as 


gtt = e 






(25) 


This is the interior time-time component of the metric tensor for the radiating 
spherical system. 

Now the tangential component of the pressure can be obtained as 


Pt = “' 55 —- Agft-m - 2g'ttgttr'^ + 2rgugttm 
d/7rr 

+2g'ttgttr'^rnr - 2gugttr^ + Agugttr'^rn + g^r^ - 2gur^m). 


(26) 


After substitution of the values of gtt , git ^nd gii , the above equation becomes 


Pt 


— \2mmr -\- Qmrmr^r — Arrn^ — 3r^mr,r 
72r^{—r + 2m)'n: 


(27) 


where mr, mr,r are the 1st and 2nd order derivative of the mass with respect to r. 


4 Physical properties of the radiation model 

4.1 Anisotropy of the system 

We have already found out the expressions for the radial as well as tangential 
pressures in Eqs. (12311 and (1271) . In terms of these pressures therefore the anisotropy 
of the system can be measured as 

— 12mmr-+ 6mrmr r — — Sr^mr-r mr 

^ = Pt-Pr = -72r2(-; + 2m)^-^ 

The nature of the pressures and anisotropy are shown in Fig. 3. which in general 
show the physically acceptable features. 
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Fig. 3 Variation of radial and transverse pressures as well as the anisotropy as a function of 
radial distance r for a = —0.1233 X 10“^ km~^ and R = 9.13 km 


4.2 Stability of the system 

4-2.1 The Tolman-Oppenheimer-Volkoff equation 

To study the stability of the system we have checked the stability equation given 
by Tolman [7] , Oppenheimer and Volkoff [8] . The TOV equation depicts the equi¬ 
librium condition of a star subject to the gravitational force, hydrostatic force and 
anisotropic force. The generalized TOV equation can be written as pm 


dr r 


where the effective gravitational mass Mg of the system is defined as 


Mg = e 2 7 . 


(29) 


(30) 


The TOV equation for our system can be translated as 

(31) 

3 gtt 3 dr r 

where the 1st term of the above equation is gravitational force, the 2nd term is 
the hydrostatic force and the last term is the anisotropic force respectively, so that 
for equilibrium of the system we should have 


Fg -|- T/j Fa — 0. 


(32) 


We have drawn the forces in Fig. 4 which describes the overall action of the 
three different forces and how the equilibrium is achieved. 


















0.004 


F. 

0.002 



Fig. 4 Variation of gravitational force, anisotropic force and hydrostatic force as a function 
of radial distance r for a = —0.1233 x 10“^ km~^ and R = 9.13 km. 

4 . 2.2 The status of the sound velocity within the system 

To examine the stability of the system here we have used the cracking tech¬ 
nique m- The condition of causality gives the physically accepted conditions 
for fluid distribution. It states that square of the tangential and radial velocities 
must lie within the limit 0 and 1. The Herrera’s cracking concept states that for 
a stable region there the radial sound speed must be greater than the transverse 
sonnd speed and also in the interior of the matter distribution two sound speed 
difference i.e. — Vsr^ should maintain it’s sign same everywhere i.e ‘no crack¬ 
ing’. In our work for the specified sets of data we find that 0 < < 1 and 

0 < Vsr^ < 1 (see Top Left panel of Fig. 5). As \vst^ — Vsr‘^\ < 1 hence our system 
is stable according to the works of Herrera m and Andreasson m- So from Fig. 
5 (rest of the two panels) it is clear that our model is perfectly stable. 


4.3 Energy conditions 

The radiating spherically symmetric system should satisfy all the energy condi¬ 
tions, viz. null energy condition (NEC), weak energy condition (WEC) and strong 
energy condition (SEC) respectively given by 


(i) NEC ■■ p + pr >0, p + pt>0, 


{a) WEC : p-hpr > 0, p > 0, p-Hpt > 0, 


(iii) SEC : p -h Pr > 0, p -|- pr + ‘2.pt > 0. 

In Fig. 6 we have shown the behaviour of all the above mentioned energy 
inequalities. 
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Fig. 5 Variation of sound speeds as a function of radial distance r for a = —0.1233 X 10 ^ 
km~^ and R = 9.13 km 



I - NEC - IVEC^ - WEC^ -sec] 


Fig. 6 Variation of different energies as a function of radial distance r for a = —0.1233 X 10 ^ 
km~^ and R = 9.13 km 


4.4 Surface Redshift 

The compactness of a star is defined as the mass-to-radius ratio of the system, 
i.e. u(r) = m(r)/r. And according to the condition of Buchdahl [13] the maximum 
allowed mass radius ratio is < 4/9 (« 0.4444) for the perfect fluid sphere. 
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For our system the compactness is 

u{r) = ^ar (35r + 252ar^ + 312a^r® + 1260ai?^ + 4368a^i?®) . (33) 

oO 

Surface red shift(^) of a star is defined as 

1 + Z= [l-2u(r)]”^, (34) 

which is for the above studied system as follow 

^ -1- (35) 

Jl - (2ar2 + + 72a?rB? + ^a^rR^) 

The compactness and surface redshift are shown in Fig. 7. From the X-ray 
spectrum the surface redshift Z can be easily observed and correspondingly com¬ 
pactness can be calculated. 




Fig. 7 Variation of compactness and redshift as a function of radial distance r for a = 
—0.1233 X 10“2 km~^ and R = 9.13 km 


5 A comparative study 

To study the physical properties of the system we choose the parameter a = 
—0.1233 X 10“^ km~^ and R = 9.13 km and mass m{R) = 1.29 solar mass (SMC 
X-1). We choose this type of data set as mentioned by Farook et al. [3] because 
our target is to investigate a radiating model which is compatible with compact 
stars. 

Using the above numerical values of different properties of the interior solution 
of spherical symmetric radiating body we have calculated each of them and also 
the energy densities of the strange star. From Fig. 6 it is clear that our model 
satisfy all the energy conditions and also other physical parametric conditions. 
We consider here a predicted star of radius 9.9 km and mass 1.71 solar mass. The 
surface density of the star is obtained as p = 6.73 x which confirms va¬ 

lidity of the model with the observational evidence mMm- However, our model 
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Table 1 Physical parameters 


Strange Stars 

Radius (in km) 

Mass (Mq) 

Mass (in km) 

PSR J1614-2230 

10.3 

1.97 

2.9057 

Vela X - 12 

9.99 

1.77 

2.6107 

PSR J1903+327 

9.82 

1.667 

2.4588 

Cen X - 3 

9.51 

1.49 

2.1977 

SMC X - 1 

9.13 

1.29 

1.9027 


cannot predict the central density as at the center there occurs singularity. The 
compactness of the predicted star is 0.2553 and corresponding redshift is 0.4295. 
This high redshift is quite relevant for strange stars. Obviously our predicted star 
have more compactness than neutron stars. Cackett et al. m have reported that 
low mass X-ray binary 4[/1820 — 30 has gravitational redshift as Z = 0.43. This 
indicates that the model studied here is a representative of a compact star and is 
suitable to explore different properties of strange stars. 


6 Conclusion 

In this paper we have tried to find out an expression of mass distribution of spher¬ 
ical symmetric radiating system by using HPM of Euler-Lagrangian equation. We 
have obtained an expression of interior solution of spherically symmetric radiating 
body and studied different physical properties of the system. 

The present investigation reveals that - 

(1) Our radiating model is compatible with the compact stars, specially that 
of strange stars as can be seen from the Comparative Study part of the previous 
Sec. 5. 

(2) From Fig. 3 we observe that the radial and tangential pressures are de¬ 
creasing with radial distance. However, up to certain distance from the center of 
the star the radial pressure becomes greater than the transverse pressure. Thus 
anisotropy is negative up to that region and after that it is increasing with radial 
distance up to the surface. This particular aspect of the last panel of Fig. 3 actu¬ 
ally shows that after 5 km, its sign changes i.e. Zi < 0 to Zi > 0 which implies that 
the anisotropic force changes from attractive to repulsive. This may be explored 
physically as follows that the fluid sphere taking its shape from oblate spheroid 
to prolate spheroid. In this connection we note reports on the non-spherical mod¬ 
els of compact stars, especially neutron stars [mfm[2oimi[2i[^ . We also note 
the recent work [23] which shows that high magnetic helds in proto-quark stars 
modify quark star masses. As we have modeled here features of the compact stars 
composed of strange quark matters so we suspect similar non-spherical effect may 
be arising in the present study also. However, it needs further investigation. 

(3) Although our model satisfies almost all the physical properties and shows 

stability of the system however seems not free from singularity. The probable origin 
of singularity may be hidden within the assumption we used in the study. In the 
total entropy of the system S of Eq. (|8|) we choose the radial component in the 
form (1 — the time component being unknown in the interior region of 

the sphere. This has been done purposely otherwise number of unknowns exceed 
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the number of equations which are not at all solvable. Since such form of the metric 
is intrinsically singnlar so there is a possibility onr metric to become infected too. 
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